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^ ■ Abstract 



C ■ We propose a procedure to extrapolate the nuclear spectral function 



P(|k|,£') obtained from nonrelativistic many-body theory to large values of 



■ three-momentum and removal energy. Our approach is based on phenomeno- 

logical information extracted from both soft hadron- nucleus interactions, in 
the regime where the proton inclusive spectrum is dictated by Regge asymp- 
totic, and deep-inelastic lepton-nucleus collisions. The extrapolated P(|k|,£^) 
is used to compute the semi-inclusive spectra of backward protons produced 

in electron-nucleus scattering. 
PACS numbers: 13.6.Le, 25.30Fj, 25.30Rw 
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Typeset using REVT^ 



The knowledge of the spectral function P(|k|, i?), giving the probability to find a nucleon 
of momentum |k| and removal energy E inside a nucleus, is a prerequisite for the theoretical 
description of a number of reactions involving nuclear targets. For both infinite nuclear mat- 
ter and light nuclei, with mass number A< 4, it is possible to carry out accurate calculations 
of -P(|k|, E) starting from a realistic nuclear hamiltonian, fitted to nucleon-nucleon scattering 
data and to the properties of few- nucleon bound states In the case of medium-heavy 

nuclei, quantitative estimates of the spectral function and its energy integral, the momentum 
distribution n(|k|), can also be obtained, using the local density approximation P]. 

The spectral functions resulting from realistic many-body calculations contain informa- 
tion on both the nuclear mean field (at low |k| and E), and short-range nucleon-nucleon 
correlations (at high |k| and E). Tipically, only about 70% of the nucleons are in the states 
of low |k| and low E, that can be described by mean field single-particle wave functions, 
while the remaining 30% are in a correlated state with another nucleon, mainly on account 
of the one-pion-exchange tensor force and the short-range repulsion of the nucleon-nucleon 
interaction. The mean field (Podkl? a^nd correlation {Pb{\^\, E)) contributions to the 
spectral function can be singled out rewriting the full -P(|k|, E) in the form 

Pi\k\,E) = Po{\k\,E) + PB{\k\,E) . (1) 

Strongly correlated nucleons play a very important role in many processes. For exam- 
ple, the production of fast backward hadrons in semi-inclusive lepton-nucleus reactions, in 
the kinematical region forbidden to scattering off a free nucleon, is mostly due to nucleon- 
nucleon correlations 0,0. While many-body calculations tipically provide a description of 
the correlation tail of -P(|k|, E) up to k 

max 

.7 GeV/c and E 

max ~ -6 GeV, the theoretical 

analysis of the available spectra of leptoproduced backward hadrons carrying large momen- 
tum P requires the knowledge of the nuclear spectral function at larger values of |k| and E 
In this note, we propose a simple phenomenological procedure that allows to extrapolate 
the correlation contribution, PB{\k\, E), beyond the region covered by nuclear many-body 
theory. 
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Let us consider a process in which a four-momentum q = (z/, q) is transferred to a nuclear 
target. Our starting point is the relationship between the nuclear spectral function and the 
function /a(-2), yielding the distribution of the nucleons in the target as a function of the 
relativistic invariant variable z, defined as 

In the above equation m is the nucleon mass, whereas Pa = {Ma,0), Ma being the target 
mass, and k = {ko, k) denote the initial nuclear and nucleon four-momentum in the target 
rest frame, respectively. The distribution function /^(-z) can be written in a general fashion 



as 



(3) 

where S{k) is the relativistic function describing the nuclear vertex with an outgoing nucleon 
of four-momentum k (see, e.g., ref. [^). S{k) can be approximated by the nonrelativistic 
spectral function according to 

S{k)=(^)p{\klE), (4) 



with 

ka = Ma- 



{MA-m + Ef + (5) 



Note that the above definition implies that fA{z) also depends upon = — g^, as pointed 
out in ip. However, in the following we will be assuming that the Bjorken limit (Q^, oo) 
is applicable, so that z can be related to the light-cone component of the nucleon four- 
momentum through z = {k~^/m) = (/cq — kz)/m {k^ = (kq)/|k||q|) and the Q^-dependence 
of /a(^) disappears. 

Substituting eq.(H) into eq.(^ we can rewrite fA{z) in terms of P{\k\,E) as 

"v^-Ma _ poo 



fA{z)=2Timz dE ci|k||k| - P(|k|,E) (6) 



where Emin is the minimum energy required to remove a nucleon from the target nucleus, 
s = {Pa + q)"^ and (for simphcity, we work in the infinite nuclear matter limit, M4 — cxd, in 
which the kinetic energy of the recoiling nucleus becomes vanishingly small) 

krain{z, E) = |m(l - z) - E\ . (7) 

The above equations imply that the large z behaviuor of /^(-z) is dictated by the high |k| 
tail of P(|k|,i?), i.e. from its correlation part Pb{\^\, E). For example, at z > 1.7 only 
momenta larger than 

kmax ~ -7 GeV/c, contribute to the integral of eq. (^). 

In ref. [Q we have developed a procedure to evaluate the asymptotic behavior of /a(^) 
without making use of the nuclear spectral function. Within this approach, based on ideas 
originally proposed in refs. [p|,p!0|, /a^z > 1) can be written as a sum of distributions fn{z/n), 
describing clusters of n (n > 2) strongly correlated nucleons, whose asymptotic behavior as 
z ^ n can be evaluated for any value of Q^. We will show that using the large z behavior 
of resulting from the approach of ref. and eq.(^ one can extract information on 

the behavior of P(|k|,£') at very large |k| and E. 

For any given |k|, the nuclear matter Pb{\^\,E) of ref. [|I| exhibits a bump located at 



-E ~ -Eyfe = y |kp + m?—m^ the energy needed to remove a nucleon of momentum k belonging 
to a strongly correlated pair of vanishing total momentum, whose width increases with |k|. 
A remarkably accurate fit to the behaviour of the correlation contribution to the calculated 



spectral function has been obtained in ref. |TT|, where Pb{\^\-iE) has been written in the 
form 

Ps(|k|,E) = n5(|k|)F(|k|,E) , (8) 

with 

POO 

nB{M)= / dE PB{\klE) , (9) 

Ethr being the minimum energy required to remove a nucleon pair. The function F(|k|,i?) 
is defined as 



4 



m{E - Ethr) - \lm{Ek - Ethr) 
T 

k 



F{\klE)=Nue^Y>{— -) , (10) 



where is related to the width through 



Tfc = 4 afc \/(2 ln2){Eu/m) = (Eb) + E^ , (11) 

{Eb) being the average removal energy associated with Pedk], E) (using the spectral func- 
tion of ref. one finds (Eb) ~ 40 MeV). The normalization constant is chosen in such 
a way as to fulfill the sum rule 

POO 

/ dE F{\klE) = l , (12) 

which in turn guarantees the overall normalization of P^dkl, E). 

Our extrapolation procedure is based on the assumption that at very large values of |k| 
and E the energy dependence of Pedkl, E) is still dominated by the contribution associated 
with the removal of a correlated nucleon pair, and can be described by -F(|k|, E) of eg. (p!0|) . 
As a consequence, the extrapolation of P{\h\,E) reduces to the extrapolation of n^dkl). 

Substitution of eq.(H) into eq.(^ leads (after inversion of the integration order) to the 
following expression for the correlation contribution to /a(-2) at large z: 

poo 

f^{z) = 2nmz / d\k\ <|.(|k|) , (13) 



with 



<|.(|k|) = |k|n^(|k|) / rfE - F(|k|,E), (14) 



implying in turn 



dz \ z 

where 

E + miz — 1 



= ^ ( fJS^ = -27rm^^{z) , (15) 



rco t -\- mi Z — 1 

^<{z)=m dE ^ — ^ ^- nB{E + m{z-l)) F{E + m{z-l),E) . (16) 

J Ethr m - E 



In general, it is not possible to invert eq.([T5|) to obtain n^dkl) in terms of (f)'{z), i.e. 
of f^i^) cind its first derivative. We have circumvented this problem using a parametrized 
ra^dkl) in eg. (p!6|) and fixing the values of the parameters in such a way as to reproduce the 
left hand side of eq.(0), evaluated according to ref. 0]. 

At kp < |k| < kjnax, kp ~ .25 GeV/c being the Fermi momentum, the behavior of the 
momentum distribution of ref. |l|] is nearly exponential, and can be accurately approximated 
using 

no(|k|) = Giexp[-(i?i|k|)-] , (17) 

with Gi = 3.30 (GeV/c)-3^ Bi = 6.2 (GeV/c)-^ and a = 1.14. 

The most natural choice is to extend the parametrization of eq.(p!7D to larger |k| and 



use it to calculate '^{z) from eq.(|l^). It should be noted, however, that the integrand in 
eq.(|T6|) has a singularity aX E = m. This problem has been taken care of inserting a cutoff 
to guarantee that n^dkl) vanish at |k| = mz. The parametrization of n^dkl) providing the 



best fit to f^iz) is 



n^dkl) = exp 



lkl - 


ko 


A- 


|k| 



no(|k|) , (18) 



where (3 = .001, ko A = mzo = m + ko + Ethr and 

no(|k|) = Giexp[-(fii|k|)"] + G2exp(-52|k|) , (19) 

with Gi = 2.90 (GeV/c)-3, Bi = 6.08 {GeV/c)-\ G2 = -21.2 (GeV/c)-^ and B2 = 26.7 
(GeV/c)~^. With the above choice of n^dkl), "^{z) is well defined and can be calculated 
using eq.(|l^) for any z > zq. 

The results of numerical calculation show that the use of eq.(|T8|) with no(|k|) given by 
eq.(|T9|) leads to a remarkably good fit to f^i^) at 2; > 1.7, the corresponding being ~ 
0.009. 

A more straightforward procedure to extract n^dkl) from eq.(|T5|) can be obtained making 
the rather drastic assumption that in eq.(|TB|) the function F(|k|,£') can be replaced by a 



5-function: 



F{E + m{z - 1), = 6{E - {yJ{E + m{z - 1))2 + - m)) . (20) 

Using eq.(iy) the E integration in eq.(0) can be readily carried out and substitution of the 
result into eq. (|l5|) leads to: 

1 (3-^^)(2-z) , 

nB(ks) = — ^-i^ d)'(z) , (21) 

2nm^ks (2 - z)2 + 1 ^ ^ ^ ' ^ ^ 

with ks = m{z — l){z — 3)/{2{z — 2)). Eq. (pl|) gives n^dkl) in terms of 4>'{z) for any values 

of z in the range 1 < z < ^/S. 

Let us now focus on the calculation of the left hand side of eq . (|T5|) . According to ref. 

@], at large values of its argument /^(2:) can be written as a sum, whose terms describe the 

contributions associated with strongly correlated n-nucleon clusters: 

/^(-) = E/n(^) ' (22) 

n=2 

the n-th term in the sum being defined for 1 < 2; < n. Within this approach the calculation 
of /^(-z) reduces to the calculation of the relevant /„(2;/n)'s, corresponding to the lowest 
values of n (tipically n = 2 and 3). The analysis of nuclear fragmentation in hadron- nucleus 
collisions carried out in refs. [^,|T3] shows that, at low Q^, the distribution of colorless three- 
quark systems in a 3n-quark cluster, Tn{z/n), exhibits true Regge asymptotic behavior as 
z ^ n. The results of refs. |9|,p!0|| provide a satisfactory description of the inclusive spec- 
tra of high-momentum protons and mesons emitted backward in proton-nucleus collisions. 
However, the sizeable Q^.^^ependence exhibited by /a (-2) at low suggests that the non- 
pertubative Q^.^lependence of Tn{z/n) has to be carefully taken into account. Starting from 
the small behavior, which can be described within the framework of Regge theory, the 
asymptotic Tn{z/n) at z ^ n and large can be obtained from the distribution of valence 
quarks inside a cluster of n strongly correlated nucleons, which can in turn be written in 
terms of the relativistic invariant phase-space volume available to a quark in a nucleon [^]. 

The function T„(z/n) can be interpreted as the distributions of effective nucleons 
within a strongly correlated cluster. Therefore, assuming that the valence quark distri- 
bution inside these effective nucleons is the same as in ordinary nucleons, the quantity 



Tn{z/n) = WnTn{z/n), Wn being the probability of finding an n-body cluster [Q, can be 
identified with fn{z/n) of eg. (p2D . The validity of this approximation, which allows to ef- 
fectively take into account nuclear excitations, is supported by the results of calculations of 
the valence quark distributions inside nucleons and baryonic resonances |[T^ , p!3[ . 



Writing the distribution of valence quarks inside a nucleon at large Q in the form 



f^^^iz) = C^z'^-il-z) 



6jv 



(23) 



where Cn is a normalization constant, qn = — ^^(O) = 1/2, a/j(0) being the intercept of the 
Regge trajectory, and ~ 2.8 — 3.2, the behavior of Tn{z/n) as z ^ n can be obtained in 



closed form |Tj]. Under the assumption Tn{z/n) = fn{z/n) we can use this result and write: 

(24) 



fn[- 



n 



n 



with = (otv + + 2)(n — 1) — 1 and A„ = a^v + ^at + 1, whereas the coefficients Dn can be 
obtained from the quark distribution in the ra-nucleon cluster, evaluated as in ref. [Q. The 
details of the derivation of eq.(p^ are given in the Appendix. 

Substituting the fn{z/nys with n = 2,3 and 4, calculated from eq. (p^) , into eq. (|2^) , one 
can easily obtain both f^{z) and its first derivative at z < 4. The resulting 0', defined as 
in eq.(|T5|), can be written: 



0'(Z) = /2 + /3 + /4 , 



(25) 



where In{z) is defined for z < n and 



I - ^ D p^^'"^ 
^ ~ 2 H2 



92-1 



- 1) 



9. I 2 



(26) 



1 /z\^3-2 

3 ^ 



1 - 



A3-2 



93-1 



^3 



(27) 



A4-2 r 



94-1 



{A, 



9. I 4 



(28) 
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Finally, the momentum distribution at large |k| can be obtained from (f)'{z) defined by 
eqs.(p5|)-p8|) using either eg. ([T5|) and the fitting procedure or the (5-function ansatz leading 
to eq.(|21|). 

Fig. |I] shows the nuclear matter momentum distribution Ti(|k|) of ref. |Q, calculated for 
|k| < .8 GeV/c (solid line), together with the values obtained from the fit to (f)'{z) at |k| > 
.8 GeV/c (diamonds). It appears that the extrapolated tail of ^(|k|) is close to a simple 
exponential dependence. The results obtained using the ^-function approximation, which 
provides an upper bound to n^dkl) at |k| > .8 GeV/ci are also shown in fig. 1 (dashed 
line). It is apparent that the energy spread of the strength is very important and cannot be 
disregarded. 

To test the extrapolated momentum distribution, we have calculated the kinetic energy 
spectrum of protons emitted at backward angle in the semi- inclusive e + A^e' + p + X 
reaction, in the kinematics of the SLAG data of ref. 0, using |^ 

PeA-.e'px{x,Q',z) = [\k\{z)] F,'' {x , Q') . (29) 

In the above equation, F<^{x,Q'^), x being the Bjorken scaling variable, is the nucleon 
structure function, while leN and leA denote the fluxes associated with scattering off an 
isolated nucleon and the nuclear target, respectively. In fig. ^ the spectrum obtained using 
the nucleon momentum distribution of ref. 0] at /c < k^ax ~ -8 GeV/c (see ref. [^) and 
setting ns(|k| > kmax) = is compared to that obtained using n(|k| > k^ax) obtained from 
the fit to /^(-z) at z > 1.7. It appears that the extrapolated n^dkl), needed to describe the 
spectrum at large proton energy (T > .3 GeV), also provides a much better fit to the data 
at lower T. 

In conclusion, we have proposed a simple phenomenological procedure that allows to ex- 
trapolate the nuclear matter spectral function beyond the region described by nonrelativistic 
many-body calculations. The main assumption needed to extract n^dkl) from the relation- 
ship between the spectral function and the distribution function /a(-z), i-e. the assumption 
that PbH^I, E), exhibits the i?-dependence associated with the removal of a nucleon be- 
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longing to a strongly correlated pair, appears to be adequate in the range of momentum 
and removal energy relevant to our analysis. The importance of the description of the in- 
dependence is clearly shown by the fact that the oversimplified 5-function ansatz leads to a 
severely overestimated momentum distribution. The numerical results shown in figs. ^ and 
suggest that our approach can be used to quantitatively investigate reactions sensitive to 
the very high momentum components of the nuclear wave function. 
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APPENDIX A: 

In ref. it has been shown that, assuming that the distribution of valence quarks inside 
an isolated nucleon can be described by 

f^S^)^z'^-{l-zf- , (Al) 

the corresponding distribution inside a strongly correlated n-nucleon cluster, obtained from 
the overlap of the phase-space volumes available to a quark inside a nucleon, takes the form 



[bN+{o-N+bN+2){n-l)] 



(A2) 



where the coefficient C„ can be written in terms of beta-functions. 

The function f^\z/n) can also be given in terms of the distribution of colorless three- 
quark systems within a Sra-quark cluster, T„(z'), using Mellin's convolution formula: 



/rm-r ft.w»fi^)#. (A3) 



where fg^{z/nz') denotes the distribution of valence quarks inside the colorless three-quark 
system. 
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Assuming that at large z/{nz') fjj^{z/nz') is approximately the same as the distribution 
inside a nucleon, e.g. assuming f^^ ~ f^, and using eq.([Al[), eg. (|A3|) can be inverted to 
obtain Tn{z') in the form: 

f^{z) = Dr,z''-{l-zy- , (A4) 

with Dn, An and gn given in terms of the constants Cn, aN and appearing in eq.( |A^) . 

Within this approximation eq. (|A^) reduces to eq.(|^), since T„ can be identifyed with 
/„, which can in turn be interpreted as the distribution of colorless three-quark systems 
inside of a strongly correlated n-nucleon cluster. 
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FIGURES 

FIG. 1. The nucleon momentum distribution in infinite nuclear matter. The sohd line cor- 
responds to the n^dkl) of ref. [Q], whereas the diamonds show the large |k| extrapolation 
parametrized according to eqs.(|l^) -(!£). The dashed line has been obtained using the (5-function 



approximation of eg. (|20D . 

FIG. 2. Kinetic energy spectra of protons emitted at angle @ in the semi-inclusive reaction 
e+CO e'+p+X. The dashed curve has been obtained using the nucleon momentum distribution 
of ref. g at |k| < k 

max ~ GeV/c and assuming Tij5(|k| > k^Y^Q^x) — Oi see ref. whereas the 
solid line shows the results obtained using the large |k| extrapolation, parametrized according to 
eqs.(p^-(19), at |k| > k^ax- The experimental data are taken from ref. p. 
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k [GeV/c] 



